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Abstract 

Properties of steady compressible flow for which geometric constraints 
have been placed on the potential function are derived, under hypotheses 
on the flow density and the singular set. Some related unconstrained 
problems are also considered, including the estimation of a class of fields 
having nonzero vorticity. 2000 MSC: 58E20, 58E99, 75N10. 



1 Introduction 

The study of certain classical fields leads to a generalization of harmonic maps 
in which the Dirichlet energy is replaced by the functional 

rQ(du) 



E= l p(s)dsdM; 
Jm Jo 



here M is a Riemannian manifold; du is the differential of a map u : M — ► N, 
where TV is another Riemannian manifold; 

Q(du) = (du,du) lT , Amu - 1TN ; 
p : M + U {0} — > M. + is a C 1,a function of Q satisfying the differential inequality 

0< p(Q) <0 ° (1) 

for Q e [0,Q cnt ). 

In the typical case, the manifold N represents a geometric constraint placed 
on the flow potential of a steady, irrotational, polytropic ideal fluid for which 
the closed 1-form du £ T*M is dual to the flow velocity. In this case we choose 

/ 7 _i \ 1/(7_1) 

p{Q)= i-^Q , (2) 
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where 7 > 1 is the adiabatic constant of the medium and Q cr it is the square of 
the sonic flow speed. These choices transform (1) into a condition for subsonic 
flow of mass density p. 

The study of hmctionals of this kind in the unconstrained Euclidean case, 
in which N = M fe and M is a domain of R", goes back at least to work on 
planar flow by Bateman in the 1920s [Ba]. An extensive bibliography covering 
the first half of the last century is given in [Be] . A more recent bibliography of 
mathematical work in compressible fluid dynamics (not necessarily connected 
with variational theory) appears in [Ch]; see also [DO] and the bibliographic 
remarks in [CF]. 

A discussion of unconstrained compressible flow in a local chart on a manifold 
appears in Sec. III. 3 of [Sed]. A global existence theorem for steady, uncon- 
strained subsonic flow on a compact Riemannian manifold is given in [SSI]; for 
subsequent research employing this nonlinear Hodge approach see, e.g., [SS2], 
[SS3], [Si], [Sm], and [ISS]. In those works the curvature of the manifold intro- 
duces geometry into the domain of the velocity field. By considering potentials 
subject to a geometric constraint, as we do here and in [Ol], we introduce geom- 
etry into the range of the velocity field. In [Ol] we emphasized this connection 
to the preceding literature by calling such potentials nonlinear Hodge maps. 
But the potentials studied in [Ol] are not associated with a cohomology class 
and neither the geometric construction nor the physical interpretation extend 
automatically to higher-degree forms. So it is perhaps more accurate to call 
mappings which are critical points of E compressible-density maps. 

There is already a considerable literature on maps for which a nonquadratic 
energy functional is given by the LP norm of the gradient of the map it; see, e.g., 
[FH], [HL], [F] and the references therein. Those works are motivated by the 
mathematical observation that the harmonic map energy is the squared L 2 norm 
of the gradient, raising the question of whether a corresponding theory can be 
derived for stationary points of the nonmctrizablc LP norms, p > 1. Our starting 
point, on the other hand, is the physical observation that harmonic maps model 
a geometric constraint on a field of constant mass density. This prompts one 
to ask whether a corresponding theory can be derived for fields having mass 
density p which depends in a nonlinear way on field strength. That question 
leads to the replacement of the harmonic map energy by the functional E. In 
applications to fluid dynamics, the variational equations of E correspond to 
continuity equations for the flow. In the incompressible (hydrodynamic) limit 
p = 1, the variational equations reduce to the harmonic map equations. This 
incompressible special case has also been interpreted in the context of nonlinear 
elasticity [T]. 

The analysis of critical points of E is necessarily somewhat different from 
analysis of the harmonic map energy or of the other LP norms of the gradient. 
In those cases scaling arguments are natural, whereas they may be unnatural 
for E, as they involve a choice of conformal behavior for p. This is a point of 
commonality between our problem and certain other recent extensions of the 
harmonic map equations, e.g., [A] and [LM]; see also [EL]. The density of maps 
which are L p -critical points of their gradient tends to zero (cavitates) as elliptic- 



2 



ity degenerates; this behavior is atypical of the mass density of fluids, for which 
the sonic value lies at the supremum of the range of subsonic speeds. [Compare 
conditions (2), above, and (23), below]. Finally, the references cited at the be- 
ginning of the preceding paragraph assume an energy minimizing property. We 
assume only that Euler-Lagrange equations are satisfied on a given subdomain. 

In Section 2 an L°° estimate is derived for nonuniformly elliptic, scalar ve- 
locity fields. Section 3 concerns technical aspects of the uniformly elliptic case 
which we studied earlier ([Ol], Theorem 3). We also present in Section 3 a 
somewhat different proof of the result in [Ol], one which is simpler in that it 
avoids certain smoothness assertions which were necessary in the original argu- 
ment. In Section 4, Corollary 8 of [Ol] is extended from the uniformly elliptic 
case to the nonuniformly elliptic case under somewhat different hypotheses. 

We note that the analysis literature tends to treat the velocity as a section 
of the cotangent bundle, whereas the physics and geometry literature puts this 
object in the tangent bundle. The local arguments of Sec. 3 are the same 
in both representations. The usefulness of the cotangent representation for 
considering fields with vorticity is apparent in Sec. 4; for consistency we employ 
this notation in Sec. 3 as well. 

2 Near-sonic maps into a line 

It is known [D] that if u e H 1 '?^) is a weak solution of the scalar equation 

div (|V?i| p ~ 2 Vu) = (3) 

for p > 1 in an open domain Q of K™ , then for every n-disc Br C f2 of radius R 
and every number S G (0, 1) there is a constant k(5) independent of R for which 

HV«||^ (BK _ 5fl) < H5)R- N I \Vufdx. (4) 

(See also [E], [Le], and references therein.) 

This is a useful result to have, as the semi-norm on the right is the energy 
integral associated to weak solutions of (3). Thus inequality (4) derives an 
L°° bound on the weak solution, which is unnatural to impose directly, from 
a condition of finite energy, which is the natural condition to impose on solu- 
tions of equations with variational structure. A uniform L°° bound on weak 
solutions plays important roles in smoothness estimates and numerical analysis. 
Inequality (4) can be significantly generalized within the class of L p -stationary 
gradients (see, e.g., [HL]). 

In this section we derive an analogue of inequality (4) for solutions of the 
scalar equation 

div [p (Q) Vu] = (5) 

for which p satisfies condition (1). Equation (5) is the Euler-Lagrange equation 
for the functional E in the special case in which u is a scalar function on R™. 
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If the middle term of (1) is bounded below away from zero on the entire range 
of values for Q, then this analogue has already been derived in considerable 
generality (see, e.g., Theorem 4.3 of [SS], Proposition 3.1 of [Si], Lemma 3 of 
[DO], or Theorem 9 of [Ol]). In each of these cases, however, the constant 
analogous to k(S) of inequality (4) tends to infinity as Q tends to Q cr it, so these 
inequalities are not uniform unless eq. (5) is itself uniformly elliptic. Rather, 
the cited inequalities contribute indirectly to uniform Holder estimates, by way 
of a delicate limiting argument introduced by Shiffman [Sh] in the planar case 
and extended to higher dimensions in [SS]. Direct arguments should suffice to 
estimate weak solutions of eq. (5) which, unlike the equations studied in the 
works cited, has scalar solutions. Our goal is to derive estimates for (5) which 
arc manifestly uniform over the entire subcritical range of values for Q. 

By a weak solution in this scalar case we mean a function it having finite 
energy E and satisfying V ip G Hq(Q) the integral identity 

/ (p(Q)du , #> dx = 0, (6) 
Jo. 

where x is a vector in a bounded type- A domain Q C K™, n > 2, and the angle 
brackets denote the euclidean inner product on 1-forms. For a definition of type- 
A domain see, e.g., p. 68 of [G]; our intention is to insure that a ball in the 
interior of Q, does not become trapped in an outward cusp. As an example, any 
Lipschitz domain is type- A In order for the following theorem to make sense in 
terms of fluid dynamics, we must additionally impose the condition that be 
topologically trivial in order that the flow potential remain single- valued. 

Theorem 1 Let the scalar function u(x), x € tt, be a weak solution in the sense 
of (6) for p satisfying condition (1). In addition, assume that p'{Q) < VQ G 
[0, Qcrit) o,nd that on this range, < k < p(Q) < p(0) < oo. Then for every 
n-disc Br strictly contained in il and every 5 G (0, 1) there exists a constant 
Ki(n,5,Ko,p(0)) for which 

sup Q(x) < Kl R- n E lBa . (7) 

x£B R(1 _ s) 

The constant n\ in Theorem 1 depends neither on the radius R nor on 
any ellipticity parameter [such as the parameter k 3 of condition (15), below]. 
Thus in particular, inequality (7) does not necessarily follow from the uniform 
bound of y/2/ (7 + 1) on the subsonic flow speed in (2). At the same time, it is 
satisfying to have a bound on weak solutions that results only from mathematical 
hypotheses on the equation itself rather than relying on a bound, such as the 
sonic speed, which is imposed on the mathematics by a physical model. The 
noncavitation hypothesis and the other hypotheses of the theorem are satisfied 
on the subsonic range by mass densities of the form (2). 

The proof of Theorem 1 is elementary The idea is to choose the test function 
in (6) to be a local restriction of the antiderivative F(Q) for the function 

f(Q)=p 2 (Q) + 2Qp(Q)p'(Q). 
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The ellipticity condition is then interpreted, where it appears, as a piece of the 
chain rule applied to the gradient of F. This relieves us of the necessity to 
bound f(Q) below away from zero, but obliges us to translate statements about 
F(Q) and its L 2 -norm into statements about Q and its energy functional. Such 
an approach combines ideas from Sec. 1 of [U] and Sec. 3 of [D]. Those papers, 
as well as [Ev] and [Le], concern weak solutions satisfying hypotheses similar 
to inequality (26) of Sec. 4. 

Proof. Wc initially assume that u is twice continuously differentiable. It is 
then easy to verify that the results are unaffected if the derivatives are replaced 
by limits of finite differences. Taking a weak derivative of (6) yields 

]T / (p(Q)u xi ) xj iP xJ dx = (8) 

where, here and below, repeated indices are summed from 1 to n. For a function 
C S Cq° {Br) and positive parameters a and (3, choose 

^(x) = u xi [Q P 2 (Q)+f3] a/2 ( 2 (x). 
Expanding the integrand of (8) yields a sum of six terms: 

n 

G°(<2K0^ = p'(Q)Qx "., "., , [Qp 2 (Q) + P] a ' 2 C 2 + 

^p'{Q)Q x iU xi u xi [Qp 2 (Q) +f3] (a - 2)/2 [p 2 (Q) + 2Qp(Q)p'(Q)] Q x3 C 

+p'{Q)Q x iu x iU xi [Q P 2 (Q) + p] a/2 2(( xJ 
+p{Q)u xixi u xixi [Q P 2 (Q) + p] a ' 2 C 2 + 
^p(Q)u xixJ u xi [Q P 2 (Q) +l3] (a - 2)/2 [p 2 (Q) + 2Qp(Q)p'(Q)] Q xJ ( 2 

6 

+p{Q)u xixi u xi [Qp 2 (Q) + p] a/2 2tt xJ =Y,*r. 

We estimate the terms of this sum individually. The following estimates should 
be interpreted as occurring "under the integral sign." The hypothesis on the 
sign of p'{Q) implies, using Kato's inequality, that 

•■ + -^(^)>([* 3 <« + ' 3 ] <o+2 " 4 )r^ 

Because the range of p(Q)/p(0) is contained in the interval (0, 1], we also have 

«2 + «5 > 
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v([Qp 2 (Q) + /3] (q+2)/4 ) 



4p(0) V^ + 2 y 
Moreover, there exists a positive constant e for which 

h + H > 

- (e P 2 /2) [p(Q) + 2Qp'(Q)f |VQ| 2 [Qp 2 (Q) + fi\ 

-2( £p 2 )- i |vci 2 [gp 2 (g) + ^] (a+2)/2 > 

2 



'c 2 . 



(«-2)/2 2 



a + 2 



v([Qp 2 (g ) + /3] (tt+2)/4 ) 



c 2 



-2(e K 2 )- 1 |VC| 2 [Qp 2 (g) + /?] (a+2)/2 . 
Choose £ to equal [k + ap(0)] /4p 2 (0). Then we obtain the integral inequality 

v([Qp 2 (Q) + /3] (a+2)/4 )| 2 C 2 *l<m / ([Qp 2 (g) + /3] (Q+2)/4 ) 2 |V C | 2 *l 

(9) 



for 



m 



2p 2 (0) 



a + 2 



«o \ K o + a/9(0) 



As a tends to either zero or infinity, m tends to a finite constant K2 which 
depends only on the upper and lower bounds on p(Q). 

Apply inequality (9.5.8) of [LU] to expression (9), taking the quantities u 
and £ of that reference to equal, respectively, the quantities \J Qp 2 (Q) + (3 and 
a/2 of expression (9). Construct a Moser iteration along the lines of expressions 
(9.5.8)-(9.5.12) in [LU]. We obtain in the limit the inequality 



xeB 



(10) 



We have 



sup [Q{x)p 2 (Q{x)) +0\< K 2 BT n f [Q P 2 {Q) + 0] * 1. 

B H(i-a) J Br 

[ Qp 2 (Q)*1= [ [ Q -^-{sp 2 (s)) ds*l< 

J Br J Br Jo aS 

[ ( Q p 2 {s)ds*l< [ p(0) [ Q p(s)ds*l = 2p(0)E lBR . (11) 
J b r Jo J Br Jo 

Regarding the left-hand side of inequality (10), condition (1) implies that Qp 2 (Q) 
is an increasing function of Q. Thus the suprema in Brm-j) of Q(x) and of 
Qp 2 {Q{x)) occur at the same value of x. Because the mass density is noncavi- 
tating, 

k 2 sup Q(x) < sup [Q(x)p 2 (Q(x))] . (12) 

i£Br(i-j) x£B r(1 _ s) 
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Comparing inequalities (10)-(12), we conclude that there is a constant K\ such 
that 

sup Q(x) +I3< K\R~ n [E\ Ba + (3vol(B R )] , 

x£B rt{1 _ 5) 

where K\ depends on n, 8, Ko, and p(0). Because (3 is an arbitrary positive 
number, we can let it tend to zero without affecting the other constants. 

We now remove the smoothness assumption on u. Replace the admissible 
test function ip(x) in eq. (6) by the admissible test function ip (x — hej) , where 
ej is the j th basis vector for E™, j = 1, n, and h is a positive constant. Then 
(6) assumes the form 

/ (p(Q(x))du(x) ,dip(x- h ej )) dx = 0. (13) 

J B 

Apply the coordinate transformation y = x — hej to eq. (13) and subtract (6) 
from (13) to obtain 

/ (p(Q(x + hej))du(x + hej) — p(Q(x))du(x) , dtp (x)} dx = 0. 
Jb 

The limiting case is an expression of the form (8). The expressions leading to 
inequality (9) remain true in the finite difference approximation. Because the 
right-hand side of (9) does not depend on h, we can allow the parameter h to 
tend to zero in this approximation, completing the proof of Theorem 1. 

Remarks. 1. In Sec. 9.5 of [LU] the Moser iteration is illustrated for 
linear equations of the form 

{Pij{x)u xi ) xi = 0, 

where 

n n 
i=l i=l 

In this case noncavitation is equivalent to ellipticity, whereas the two conditions 
are distinct for the quasilinear density p(Q). Thus the ratio p/v in expression 
(9.5.8) of [LU], which is analogous to the factor m in our expression (9), intro- 
duces a dependence on ellipticity in the linear case but not in the quasilinear 
case. 

2. Theorem 9 of [Ol] is a subparabolic Moser estimate for multivalued flow 
potentials possessing geometric constraints. The preceding proof is too simple 
to work in that case, and the constants obtained in the proof of Theorem 9 
depend on ellipticity. However, one can replace, in the line preceding inequality 
(69) of that proof, the function Q r_1 for r > 2 by the function (Q + /3) r_1 for 
r > 1, allowing (i to tend to zero at the end as in the preceding proof. This 
avoids eventual difficulties in the Moser iteration. 
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3 Uniformly elliptic maps 



We now consider the more difficult cases in which the target has nontrivial 
geometry. In what follows the symbol C will denote generic positive constants 
unless otherwise indicated. 



3.1 Effects of geometric constraints 

In studying critical points of E, it is natural to obtain the admissible class of 
maps from the condition of finite energy. We seek a class of bounded maps 
having integrable density e{u). In order to integrate this object, it is necessary 
to choose local coordinates for e(u) on N and it is not a priori true that u takes 
a coordinate chart on M into a coordinate chart on JV. If however we restrict 
our attention to maps from M into N which are Holder continuous, then the 
local oscillations of the map are controlled on the target, and the image of a 
sufficiently small region of M will lie in a coordinate chart of TV. In this case we 
can write 

where for n = dim(M), x — (x 1 , x n ) is a coordinate chart on the manifold M 
having metric tensor -f a/3 (x)- 1 u = (u 1 ,...,u m ) is a coordinate chart on the 
manifold N having metric tensor g i j{u) 1 where m = dim(N); repeated Greek 
indices are summed from 1 to n; repeated Latin indices are summed from 1 to 
m. 

This continuity assumption severely restricts the kinds of questions that we 
can ask about the map. Moreover, the geometric constraint re-emerges as a 
problem, even if the map is continuous, when we attempt to extremize the 
energy functional by taking variations. This is because the test functions ip 
might take the image of u + tip off of N, even for small values of t. This can be 
immediately seen if, for example, we take N to be the unit sphere |u| = 1. 

The conventional solution to both problems, that of defining an admissible 
class of finite-energy maps and of varying the energy on the target manifold, 
is to embed N isometrically into a higher-dimensional Euclidean space R fc by 
the Nash Embedding Theorem. The manifold N emerges as a system of k — m 
independent constraints, 

$(«) - (*l(ti),. :,$k-m («)) =0. 

In this case 

and the incompressible energy integral reduces to the classical Dirichlct integral. 
In taking variations, a suitable Euclidean neighborhood 0{N) of N is projected 
onto N by nearest point projection II. If t is small enough and N is a C 1 
submanifold of M fc , then the variations ITo (u + tip) will be constrained to lie on 
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N for almost every x in M, where ip £ C§° (M, R fe ) , and for every x in M if u 
is continuous. Now the variational equations of E are given by 



dt\t=oJ M J 




where 



Q t = \d[Uo{u + tiP)]\ 2 . 



The variational equations in the ambient space assume the explicit form 



where 5 is the formal adjoint of the exterior derivative d and A is the second 
fundamental form of N. 

See [Sch] and [ScU] for detailed discussions of these issues in the harmonic 
map case. 

We call u e L°° (M,R k ) n H 1 * 2 (M,R k ) a weak solution of eqs. (14) in a 
coordinate chart fl of M if u has finite energy E and satisfies, V < G #o' 2 (fi, Rfc )n 



The existence of weak solutions to the unconstrained problem in the elliptic 
range follows, by lower semicontinuity, from the convexity of the energy func- 
tional under condition (1). Weak solutions of the constrained problem may not 
exist for certain choices of p and N. To see this, let p(Q) = Q^ p ~ 2 ^ 2 , p > 1, 
and consider the counterexample of [HL], Sec. 6.3. 

3.2 Maps with apparent singularities 

The literature on removable singularities is too large for even a superficial re- 
view. We mention that the removability of singularities in harmonic maps is 
considered in, e.g., [SaU], [EP], [Me], [Li], and [CL]. Obstacles to the exten- 
sion of methods used in those references to our case include, in addition to the 
dependence of the scaling behavior of E on the choice of p, the absence of an 
obvious analogue to the a priori Holder estimate of [HJW], which forms the 
basis for many smoothness results in the harmonic map literature. Removable 
singularities theorems and related a priori estimates for mappings which are 
critical points for the L p -norm of their gradient are reviewed in [F] . Those argu- 
ments also strongly depend on the scaling behavior of the energy. Removability 
of singularities in systems which resemble the unconstrained case of eqs. (14) 
can be found in, for example, [ISS]. The application of such results to the con- 
strained case is limited by the presence of quadratic nonlinearities arising from 
the target curvature. 

The removability of an apparent singularity can be proven either by showing 
the existence of a continuous transformation to a nonsingular domain, or by 



S [p (Q) du] = p{Q)A (du, du) , 



(14) 



L°° (n,R k ) , the identity 
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ruling out the existence of the singular set on a priori grounds. We adopt the 
latter approach in the following theorem. 

Theorem 2 Let u : Q — > N be a C 2 stationary point of the energy E on ri/£, 
where is an open bounded, type- A domain of W 1 , n > 2; N is a smooth, 
compact m- dimensional Riemannian manifold, m < n, dN = 0;EccBccf! 
is a compact singular set, completely contained in a sufficiently small n-disc B, 
which is itself completely contained in tt. Suppose that p satisfies 

k 3 < p(Q) + 2Qp'(Q) < k 4 (15) 

for constants < k 3 < m < oo. Ifn > 4, let 2n/ (n — 2)<p<n, where p is the 
codimension o/S, and let du G L n (B); ifn = 3,4, let du G L 4 i 0,3 (B) n L 4 i(B), 
where (3 = (p — e) / (p — 2 — e) for 2 < p < n, e > 0, and | < q < q. Then du 
is Holder continuous on B. 

Because the singular set is assumed small, the choice of a Euclidean domain 
entails little reduction in generality. In distinction to the harmonic map case, 
Theorem 2 does not immediately imply any higher degree of smoothness. The 
theorem immediately extends to the case of a finite number of singular sets 
having the same properties as X. 

Theorem 2 is stated and proven in [Ol] (Theorem 3). We begin by briefly 
reviewing that proof, adding details on the underlying elliptic theory in Lemmas 
4 and 5. An alternate method of proof, which avoids Lemmas 4 and 5 altogether, 
is given in Sec. 3.3.1. We show in the proof that the modulus of continuity for 
du depends on p, u, N, n, and on the I/ n -norm of du. A metric can be chosen on 
in which the L"-norm of du over f2 is smaller than any given fixed number. 
There are choices of p, however, under which the variational equations fail to 
be invariant under this transformation; c.f. [KFL]. 

Lemma 3 Under the hypotheses of Theorem 2, u is Holder continuous on B. 

Proof. Away from the singular set, Q is sufficiently smooth that local coor- 
dinates can be chosen on N and one can show ([Ol], Theorem 2) that 

L(Q) + C N Q (Q + 1) >0, (16) 

where L is an elliptic operator under hypothesis (15). Integrate inequality (16) 
by parts over B against a test function (ntjj) S(Q); here rj, > 0; ip(x) = V.x in 
a neighborhood of G C§°(B') where E CC B' CC B; E(Q) = h{Q)ti (Q) , 
where for k = 0, 1, . . . , 

h(Q) = 

Q[n/(n-2)] k n/i f O r0<Q <£, 



forQ> 



if n > 4; h is an analogous test function ([Se], p. 280) when n is 3 or 4. Let tp be 
the limit of a sequence 1 — where £^ is the sequence rf v>> of [Se], Lemma 
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8. This sequence has the property that £^ — a.e. in a neighborhood of S, 
but £^ tends to 1 a.e. and tends to zero in L^~ 6 as v tends to infinity. 

It can be shownQOl], (28)-(35)) that these choices imply the inequality 

J v 2 \v(Q T(k)/4 )\ 2 *l< J Vr7| 2 Q" (fc)/2 * l, 

where r(fc) = n[n/(n — 2)] k . Taking k to equal zero, the right-hand side of 
this expression is bounded by the L n -norm of du over Q. Applying the Sobolev 
inequality to the left-hand side allows us to repeat the preceding integration 
by parts for k = 1. Applying the Sobolev inequality to the resulting inequality 
allows iterations for progressively higher values of k. In this way any finite L p - 
norm for du can be bounded by the L"-norm of du over il. We conclude that 
du lies in the space L P {B) for any finite value of p and is an H 1 ' 2 weak elliptic 
subsolution on B' . Then u is Holder continuous and the proof is complete. 

Let D be an n-disc of radius i?, completely contained in the n-disc B', 
completely containing the singular set £ and centered at a point io G S. We 
require a classical result on linear boundary- value problems: 

Lemma 4 If p € C 1 ^ (D) and w £ C° (3D) , then 3v : D -> R fe such that 
v G C 2 ' a (D) n C° (TP) and v satisfies the linear boundary-value problem 

p(\x\ 2 )dv\ = inD R (x ), (17) 
v$ = Wtf on dD, 

where the subscripted ■& denotes the tangential component of the map in coordi- 
nates (r, #i, ...i? n _i). 



Proof. Condition (15) implies [U] that 

> K 



p(M 2 ) 



for some positive constant K. This inequality implies strict cllipticity of the 
linearized equations (17). The result now follows from Theorem 6.13 of [GT], 
although that result is stated for scalar equations, because the differential op- 
erator in (17) is diagonal. This completes the proof. 

Define a map <p : D — > R k and consider the nonlinear boundary-value prob- 
lem 

5 (p(\ckp\ 2 )chp^ = inD R (x ), (18) 
= u$ on dD. 

Lemma 5 If u satisfies the hypotheses of Theorem 2, then the boundary-value 
problem (18) has a solution in the space C 2,/3 (D) fl C°' a (TP) . 
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Remarks. If u lies in the space C 2 '' 3 (D/H) and if £ is an isolated point (or 
by extension, a finite point set), then the smoothness of p follows from Schauder 
estimates, and as the radius R of D shrinks to a point, the boundary conditions 
of problem (18) remain smooth. In this case one can compare du to dip with 
the goal of applying Theorem III. 1.3 of [G] exactly as in [Ol], and no further 
remarks are necessary. If, however, £ is not a point, then for sufficiently small 
R, D r will intersect S and we have only the result of Lemma 3, that u is Holder 
continuous on a domain that includes the singular set. It is not explicitly shown 
in [Ol] that this is sufficient boundary regularity for completing the proof of 
Theorem 2; but that is in fact the case, as we will show here. 

Proof of Lemma 5. Consider the boundary- value problem (17) for iu = Mon 
D r (xo), where r £ (0, R\. The boundary data are Holder continuous by Lemma 
3, so the solution v lies in the space C 2 ' a (D r ) n C° (D r ) by Lemma 4. Now 
we extend to systems the proof of [LU], Theorem 4.8.7. That is, we solve a 
sequence of boundary- value problems having the form 

5 [p{\dp\ 2 )dp^j = in D rt (x ), D ri C D R , (19) 
<p# = v$ on dD ri , 

where {r{\ — > R. A C 2,a solution p i to this problem exists for every i by 
Theorem 1 of [SS2]. (The differentiability requirements on the boundary are 
encapsulated in the definition of the space T>2 of that paper.) Also, by hypothesis 
u is bounded by a constant depending only on N. This gives a uniform bound 
on the boundary data vs on each dD ri . Solutions <p i of (19) satisfy a maximum 
principle, for each i, by Sec. 2 of [SS3]. Thus the sequence {p { } possesses a 
subsequence which converges, as j-j tends to i?, to a solution ip € C 2 ^ (Dr) n 
C°- a (Dr) , as required. This completes the proof of Lemma 5. 

We now complete the proof of Theorem 2 by showing that the differential 
du is Holder continuous in a domain that includes the singular set. 

For sufficiently small B, we can construct a suitable n-disc, on the boundary 
of which the tangential component of a comparison vector can be forced to agree 
with the tangential component of u (c.f. [Li], Sec. 3). 

Consider a solution p to the boundary-value problem (19). Combining 
Lemma 5 with Theorem III. 1.2 of [G], we find that if (dp)R. X() denotes the 
mean value of the 1-form dp on Dr(xq), then for any sufficiently small R, dp 
satisfies 

f \dp-(dp) R , XQ \ 2 *l<CR n+2X 

JD R (xo) 

for some number A G (0, 1]. Then 



Dr(x ) 



d(u — p), p(\du\ 2 )du — p(\dp>\ 2 )dip 



* 1 



Dr(xo) 



(u — p) , p(\du\ 2 )A(du, du) ) * 1. 
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Apply a generalized mean-value formula to the 1-form p(\ds\ 2 )ds as in [Si], 
Lemma 1.1. We obtain 

/ \d(u- (p)\ 2 *1 < 

JD R (x a ) 

{\du\ + \dip\) \x\ * 1 + f \u-ip\ p{\du\ 2 ) \u\ \du\ 2 * 1 ) 



lD R (x ) 

= ii+i 2 - (20) 

We have 

h 



< [ {\d(u-ip)\ + 2\ckp\)\x\*l< 

JDr(x ) 

e \d(f- u)\ 2 * 1+ 

JD R {x ) 

R 

n+1 



x\ n+L d\x\, 



where the sup norm of dip depends on the modulus of continuity for u through 
cq. (18). 



h < C (p) \ \u — tp\ \u\ \du\ 2 * 1 < 

Jd r (x ) 

R~ v [ \u - ip\ 2 \u\ 2 * 1 + R v [ |du| 4 *l< 
R~ v I \u- V \ 2 \u\ 2 *l + C (\\du\\l , n) R"(p-i)/p+» 

JDr(x ) V ' 



'Dr(xo) 

for a constant v to be chosen and p so large that vp > n. We have by the 
Sobolcv Theorem 



R~ v [ \u - ip\ 2 \u\ 2 * 1 < 

JDr{x ) 

/ . \ (n-2)/n , x 1 , 

/r / | u -^| 2 " /( - 2) *i / h"*i 

\Jd r (x ) J \JDr(xo) 

<R- V C S obolev I \d{u-p)\ 2 *l ( [ 

JDr(x ) \Jl 

[ \d(u-^)\ 2 *l[[ 
Jdr(x ) yo 



\ 2/™ 

i«r*i ) 

\ 2/n 



<R- v C\\u 2 \\„ n „„^ / \d(u-w)\ 2 *l / \x\ n - l d\x\ 



<R 2 - U C / \d{u-ip)\ 2 *l. 
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Choose v e (0, 2) . Substitute the estimates for i\ and 12 into the right-hand 
side of (20) and absorb small terms on the left in (20) to obtain 

/ \d (u - if) | 2 *1 < CR n+y 

Jd r {x ) 

for some positive A'. 

The minimizing property of the mean value with respect to location param- 
eters implies 




< / \du-dip\ 2 *l + I \dcp - (dip) R . Xa \ 2 * 1 

Jd r (x ) Jd r (x ) 

< CR n+e (21) 

for some £ > 0. Because these estimates can be repeated for any sufficiently 
small value of R, the proof of Theorem 2 is completed by the local form of 
Campanato's Theorem (Theorem III. 1.3 of [G]). 

3.3 Weak solutions of eqs. (11) and (15) 

If we make no assumptions about the singular set but assume that u satisfies 
(14) weakly in B, it is possible to show by estimating difference quotients that 
du is an element of the space H 1,2 (B). The next logical step would be to show 
du € L°° (B) . This step cannot be taken in the constrained case by following 
an analogy to the unconstrained case. The latter arguments proceed from a 
scalar inequality, as in Sec. 1 of [U]; but in order to adapt that argument to 
the constrained case it is necessary to choose local coordinates on the tangent 
space of N, as in Theorem 2 of [Ol]. This requires some a priori information 
about the singular set of u. 

Theorem 2 of the preceding section implies that if an _B-critical map u is 
bounded and Holder continuous on an open Euclidean domain, then du is Holder 
continuous on small compact subdomains. The initial continuity assumption 
means that the target geometry will play little role in the analysis beyond its 
contribution to the nonlincarity of the variational equations. 

3.3.1 An alternate proof of Theorem 2 

The arguments of [U] imply that weak solutions of (18) are Holder continu- 
ous, but the Holder estimate implied by that work cannot be continued up to 
the boundary. Nevertheless, it is possible to show that interior smoothness of 
weak solutions to (18) is sufficient to complete the proof of Theorem 2 using a 
modification of the preceding arguments. 

The weak form of problem (18) can be written 
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L 



dC,p(|#| 2 ) dip) * 1 = 0, (22) 



'-Dr(xo) 

where d( is a closed 1-form in L 2 (D R ) having vanishing tangential component 
on dD R . Applying standard function-theoretic arguments on M™, we consider 
( to be an admissible test function; c.f. eq. (1.2) of [Si]. Writing (22) as the 
weak variational equations of the energy functional E with TV replaced by R fe , 
we have 

f f\ d v\ 2 f , 

/ / p(s)ds*l>K \dip\ 2 * 1, 

JD R (x )Jo Jd r (x ) 

so dip lies in the space L 2 {B R ) by ellipticity and finite energy. (See also Sec. 1 
of [U].) The proof of Lemma 3 implies that du lies in L 2 (D R ). Because d(u — ip) 
is in L 2 , we can choose ( = u — ip in (22). The resulting weak Dirichlct problem 
is solvable by Proposition 4.3 of [Si]; see also [ISS]. The 1-form dip is Holder 
continuous in the interior of D by Proposition 4.4 of [Si], which is derived from 
[U]. The Campanato Theorem implies that 

f \dv-{d v ) R/2tXQ \ 2 *l<CR n+a 

for some a £ (0, 2]. Estimating (20) as in the preceding section, we find that 

\d(u-p)\ 2 *l <CR n+f * 



L 



'D R (x ) 

for some positive /x. Then of course 



f \d(u-ip)\ 2 *l<CR n+ ». 

J D R/2 (x ) 

Rewrite inequality (21) over D R / 2 (xo) to obtain 

/ \du - (du) R/2yXa \ 2 * 1 < / \du - {dip) R/2 , Xo \ * 1 

Jd r/2 (x ) Jd r/2 (x ) 

< \du-dip\ 2 *l+ \dip - {dip) R/2 , XQ \ 2 * 1 

Jd r/2 (x ) JD a/2 (x ) 

< CR n+e 

for some I > 0. This completes the alternate proof of Theorem 2. 

An application of this argument to an unconstrained problem for bundle 
curvature is given in [02]. 
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4 Nonuniformly elliptic solutions having nonzero 
vorticity 

Note that the map u enters into the problem of the preceding section only 
through its geometry. The variational equations in the unconstrained case are 
written in terms of du, and the map u does not directly appear in them. This 
raises the question of how much of the theory can be deduced in the uncon- 
strained case without assuming the existence of a potential. 
Thus we consider systems having the form [Ol] 

5(p(Q)w) = 0, (23) 

duj = v A uj, (24) 

for v € A 1 (V A ), where V is a smooth section of a vector bundle over an open, 
bounded domain Q of R"; uj e A p (V) ; Q = uj A *uj, where * : A? A™"? is the 
Hodge involution; p is defined as in the preceding sections, but will be assumed 
to satisfy an inequality somewhat different from (1). 
The condition 

duj = (25) 

implies, by the converse of the Poincare Lemma, the local existence of a potential 
u G A p_1 (V) such that du = uj. Thus solutions of (25) lie in a cohomology 
class, whereas solutions of (24) do not in general. However, the integrability 
condition (24) generates a closed ideal when p = 1. Obvious modifications of 
condition (24) generate a closed ideal for solutions of higher degree (see, e.g., 
[Ed], Theorem 4-2.1). 

In the interpretation of eqs. (23), (25) in which w is dual to the flow velocity, 
the vanishing of the vorticity dui expresses the property that uj must integrate 
to zero along any curve homologous to zero. The 0-form u is the flow potential. 
One-form solutions to (24) have only the weaker property that uj = idu for some 
0-form t. 

Remark on terminology. The distinction between curl-free and rotation- 
free fields is sometimes used to characterize velocity fields corresponding to (25) 
and (24), respectively {c.f. pp. 123, 124 of [MTW] ; p. 28 of [So]). In [Ol] 
the term irrotational field is used to denote a curl-free field, and any field which 
is not curl-free is called rotational. While that terminology may be misleading 
physically, the term curl-free is not mathematically correct in higher dimensions, 
so cither choice of terms is open to criticism. 

In Theorem 7 and Corollary 8 of [Ol] a Holder estimate is derived for the 
variant (23), (24) of the nonlinear Hodge equations on a possibly singular do- 
main. As the solution approaches the critical value at which the ellipticity of 
the differential operator breaks down, the elliptic estimate of [Ol] also breaks 
down. In this section we derive an estimate which is uniform over the entire 
subcritical range. 



16 



We assume that w is a classical solution of eqs. (23), (24) outside a singular 
set of prescribed dimension and that the density p satisfies 

+ k Q ) q < p(Q) + 2Qp'(Q) < k 5 (Q + k ) q , (26) 

for constants K$,q > and fco > 0. Condition (26) was introduced in [U] in 
connection with a generalized version of eqs. (23), (25). That condition implies 
that there is a possibly larger value of K5 for which 

^\Q + k y <p(Q) <K 5 (Q + fco) 9 (27) 

and 

\Qp'{Q)\ <K 5 (Q + k)". (28) 

In the sequel we denote by n a number so large that it satisfies (26), (27), and 
(28). Condition (26) is an ellipticity condition for eqs. (23). If k — 0, then 
ellipticity degenerates as Q tends to zero; condition (27) implies that the density 
p also tends to zero (cavitates) in this limit. Thus ellipticity and noncavitation 
are equivalent under condition (26). In applications to compressible flow, the 
degeneration of ellipticity need not imply cavitation, and in cases in which 
these two phenomena are equivalent, as in the Chaplygin approximation, the 
degeneracy occurs at infinity rather than at zero. Moreover, condition (26) is 
not associated with a sonic transition. For these reasons, condition (26) does 
not appear to be appropriate for applications to fluid dynamics. However, it 
arises in certain natural generalizations of the Dirichlct energy; see [HL] and 
the references cited therein for details. 

The methods used to study eq. (23) also apply to systems in which (23) is 
replaced by an equation of the form 

6(p(Q)u)=t(u>), (29) 

where p = 1 and £ is a scalar function of u satisfying 

\?(u>)\<K(Q + k) a (30) 

for a E K+. For simplicity we take a = q. Obvious algebraic modifications lead 
to results analogous to inequality (34), Theorem 6, and Corollary 7 for general 
a > 0. In that case, inequality (34) may no longer be linear in its terms of zero 
order. 

Certain properties of eqs. (29), (24) can be obtained by deriving a differential 
inequality for an appropriate scalar function of the solution. The case v = £ = 
is due to Uhlenbeck [U], who framed the argument in the context of a broadly 
defined elliptic complex. We initially present a version of Uhlenbeck's argument 
in simpler notation for solutions of (29) and (25), and then indicate how to 
modify the proof for the case of solutions of the system (29), (24). 

Denote by H(Q) a C 1 function of Q such that 

H'{Q) = \p(Q) + Qp'(Q). 
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Then [U] 

(u, A (p(Q)w)) = (u, d t (p(Q)uj)) - (diuj, di (p(Q)w)) 

= AH(Q) - \p{Q) (d t io, duo) + p'(Q) w) 9,Q] , 
where di = d/dx % , x = x 1 , x n <G ft, and 

p'(Q) (diW,w)diQ = J2 2 P'(Q) {^,lo) 2 . 

i 

We have [for either sign of p' (Q)] 

p(Q) (doj, diu) + p\Q) (8 1 lo, lo) d l Q > k-\Q + kf \Vu\ 2 

and 

(w, A (p(Q)w)) < Air(Q) - k-^Q + |Vw| 2 . 

In addition, 

(w, A (p(Q)w)) - (w, «Jd (p (Q) co)) + (w, ctf (p (Q) w)) 

= (u,6d(p(Q)uj)} + (uj,dO 
for solutions of (29) and (25), yielding 

k- 1 (Q + k) q \Vuo\ 2 < AH(Q) - {lo, Sd (p (Q) lo)) - (lo, d£) . (31) 

Define a map /3 W : A — > A p+1 by the formula /3 U : ( —> d( A lo, for ( E A 
and lo E A p . Ylv = 0, 

U (C) = d(Cw), 

but we do not use this property. Define the map /3* : A p+1 — > A by the 
formula 

/?* (?) = * * (w A *?) 

for?eA p+1 . Writing 

(w, <M (p (Q) w)) e*[wA *(5d (p (Q) lo)} 

= *d[wA* (p' (Q) (IQAw)]=(5*[wA* (p' (Q) rfQ A w)] 
(c.f. Lemma 2.1.4 of [J]), we can write (31) in the form 

k~\Q + kf |Vc| 2 < AH (Q) - flX [p (Q)] - ( W) d0 . 

Because 

we can rewrite this inequality, in terms of H, as 
k _1 (Q + fc) 9 |Vw| 2 < AH - di« |* 
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w A * 



P 7 (g) 

(<3) 



dH Alo 



for 



Write 



AH -ft [q u (H)]-{u>,d£). 



(H) = AH -ft, [?„ (H)} 

k,j 

If p' (Q) is nonpositive, then the matrix dkj satisfies 

i<„ i i W(g)j <i i MQ + fc) 9 



ff'(Q) " k- 1 (Q + k)i' 
Letting ir = (ni, . . . , ir n ) denote an n- vector, we have 



k| 2 < ^2 KkakjKj < (l + 2k 2 ) |tt| 2 . 



If // (Q) > 0, write 



div 



1 



Op' (Q) 



&(Q) + Qf/(Q)-Qf/ (Q) 



dw 



(Q) 
p(Q) 



grad (H) 



2H> (Q) 



grad (id) 



The matrix a,kj now satisfies 

- 1 p(Q) 



< 



2H'(Q) 

— Clkj < -r^TTTTTT = 1. 



2k ~ 2H> (Q) " 2id' (0) 

Letting tt = (m, . . . , ir n ) denote an n- vector, we have 

I 71 "! 2 ^ i i 2 



Thus _L is a uniformly elliptic operator on id for either sign of p' (Q) . 

It remains to estimate the lower-order nonlinear term (u>, d£) and to adjust 
for v ^ 0. We have 



(w, d£) = (u, £' (lu) duj) = £' (w) (uj, v A oj) 



> 
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- 1£' M| W\Q>- |£' H| W\ (Q + k)> -k \v\ (Q + k) q+1 , 

using (30) with a = q. Integrating condition (26) over a dummy variable in 
[0, Q] and using H(0) = 0, we find that 

K- 1 {Q + ky +1 <H(Q) <K{Q + kf +1 (32) 

and obtain 

(w,d£) > -k 2 \v\ H(Q). 

In the proof of Theorem 7 of [Ol] it was shown that if lo, v smoothly satisfies 
(23), (24), then there is an independent positive constant C and a sufficiently 
small constant e (k) for which 

< (k- 1 - ek)(Q + k) q \Woo\ 2 

< AH(Q) - *d [lo A * (p'(Q)dQ A u)] 

+ C(Q + ky (\Vv\ + \v\ 2 )Q. (33) 
We can convert this estimate to an inequality in H, noticing first that 

(0 + k)« (\Vv\ + \v\ 2 ) Q < k (\Vv\ + \v\ 2 ) H (Q) 

by (32). Now taking into account the term — {us,dS,) and reasoning as in the 
curl-free case, we rewrite (33) in the form 



< L u (H) + C (k, q) (|Vv| + \v\ 2 + t\v\) H, (34) 



where 



L w (H) =AH- div 



■ 1 WW) dHA ' 



and t = unless £ is nonzero, in which case t = 1. This operator is clearly 
elliptic on _ff , as we did not use the closure of lo under d in establishing uniform 
cllipticity for the corresponding operator in the case v = 0. 

Notice that the operator L w (H) can be written as an operator on Q having 
the form 

L u {Q) = di {^p(Q) + Qp l (Q^jd i Q 

-*d[w A*(p'(Q)dQ Aw)]. 

This operator is only elliptic on Q only if k exceeds zero. Thus for example, 
inequality (34) allows us to extend Corollary 8 of [Ol], which was based on 
an elliptic inequality for L u (Q). The bound on lo established in that result is 
not uniform as Q tends to zero unless the constant k in condition (26) exceeds 
zero. We can remove that restriction if we place different LP hypotheses on the 
solution. In comparison with the hypotheses of [Ol], Corollary 8, the new LP 
hypotheses placed directly on lo are somewhat stronger, whereas those placed 
indirectly on lo, through the hypothesis on v and its derivatives, are considerably 
weaker. 
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Theorem 6 Let the pair to, v smoothly satisfy eqs. (23), (24), with p satsifying 
condition (26), on a domain Q/E. Let CI be a type- A domain ofW 1 , n > 2. Let E 
be a compact singular set, completely contained in a sufficiently small n-disc B, 
which is itself completely contained in CI. If n exceeds 4, let 2nj (n — 2) < /! < n, 
where \x is the codimension of E, and let to lie in L™(« +1 ) (B) . If n — 3 or 4, 
let u lie in L 4 («+ 1 )/ ? 7i ( B ) n L 4 («+ 1 )T2 (#) , where j3 = (p — e) / (fi — 2 — e) for 

2 < jU < n, e > 0, and ^ < -f 1 < 7 2 . If \Vv\ + \v\ 2 6 L s (B) for some s exceeding 
n/2, then lo is bounded on compact subdomains of CI. 

Proof. Integrate inequality (34) against the Serrin test function as in Lemma 

3 of the preceding section. Using (32), the L p hypothesis on lo translates into 
LP hypotheses on H which are sufficient for applying the arguments of Lemma 
3 to H. These yield an integral inequality which can be iterated. After a finite 
number of iterations, we find that H is in L p for all finite values of P and is 
a weak H 12 subsolution on B n E. Theorem 5.3.1 of [Mo] implies that H is 
bounded on compact subdomains of CI. Condition (26) extends this result to 
Q, and thus to lo. This completes the proof of Theorem 6. 

Corollary 7 Assume the conditions of Theorem 6, except let lo be a 1-form, 
replace eq. (23) with eq. (29), and let £ satisfy (30). Then the conclusion of 
Theorem 6 remains valid. 

Proof. Clearly, |Vu| + \v\ 2 + \v\ e L s (B) for some s exceeding n/2. This 
completes the proof of Corollary 7. 

Remark. We take this opportunity to correct a pair of misprints in the 
statement of Corollary 8 of [Ol]: replace L 4q (B) by L 4qi (B) and 1/2 < q < q 
by 1/2 < q < qi. 
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